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ABSTRACT 

The central engine of gamma-ray bursts (GRBs) is believed to eject double-sided ultra-relativistic 
jets. For an observed GRB, one of the twin jets should point toward us, and is responsible for the 
prompt gamma-ray and subsequent afterglow emission. We consider in this Letter the other receding 
jet, which will give rise to late-time radio re-brightening (RRB) when it becomes non-relativistic (NR) 
and radiative isotropic. The RRB peaks at a time 5tNR = 2{Ej^^i/nY^^ yr after the GRB, where tNR 
is the observed NR timescale for the preceding jet, Ej is the jet energy and n is the ambient medium 
density. The peak flux is comparable to the preceding-jet emission at Inr- We expect the RRB of GRB 
030329 1.7 yr after the burst with a flux ~ 0.6 mJy at 15 GHz. The cases of GRBs 970508 and 980703 
have also been discussed. The detection of RRB, which needs dense monitoring campaign even a few 
years after a GRB, will be the direct evidence for the existence of double-sided jets in GRBs, and prove 
the black hole-disk system formation in the cores of progenitors. 

Subject headings: gamma-rays: bursts — ISM: jets and outflows — radio continuum; general 



1. INTRODUCTION 

A standard shock model of gamma-ray burst (GRB) 
afterglows has now been well established, in which the 
GRB outflow (jet) drives a relativistic expanding blast- 
wave that sweeps up and heats the GRB ambient medium 
to produce the long-term X-ray, optical and radio after- 
glow by synchrotron/inverse-Compton (IC) emission (see 
Zhang & Meszaros 2004 or Piran 2004 for recent reviews) . 

The remarkable supernova signature in GRB 030329 af- 
terglow (Hjorth et al. 2003) has confirmed that GRBs, 
at least the long-duration ones, are originated from ex- 
plosions of massive stars. Their progenitors are be- 
lieved to form black hole-disk systems in their cores, and 
produce double-sided jets, which are responsible for the 
GRB emission, along the spinning axis of the black hole 
(Meszaros 2002). For the observed hard GRB photons 
to escape freely, avoiding electron/positron pair forma- 
tion, the outflows are reguired to be ultra-relativistic with 
Lorentz factor F > 100. Because of the relativistic beam- 
ing effect, the GRB emission is confined in a narrow cone, 
and can only be observed when the light of sight is within 
the cone. Interestingly, it was this beaming effect that 
made many authors expect the present of "orphan after- 
glows" (Rhoads 1997; Dalai, Griest, & Pruet 2002; Gra- 
not et al. 2002; Granot & Loeb 2003) which correspond to 
off-axis GRBs without gamma-rays associated, but with 
afterglows detected when the shocks decelerate. 

The previous works, except for Granot & Loeb (2003), 
only discussed the emission from the preceding jet (PJ), 
which is pointing toward us, and neglected that from the 
other receding jet (RJ), since it always points away from 
us. However, in this Letter we discuss specially the emis- 
sion from the RJ, which will emerge when the RJ become 
non-relativistic (NR) and the radiation becomes isotropic. 
Because of the light travel time delay, the RJ contribu- 
tion will overlay above the late-time, already decayed PJ 
emission. In the NR phase the afterglow peak emission 
has moved to radio frequencies, so the re-brightening from 
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the RJ is only expected in late-time radio afterglows. We 
predict the late-time radio re-brightening (RRB) of GRB 
030329 and discuss the long-term radio observations for 
two GRBs 970508 and 980703. The observation of RRBs 
will be particularly important due to its straightforward 
implication of the existence of double-sided jets in GRBs. 

2. EMISSION FROM THE RECEDING JET 

Consider two collimated jets, that are ejected in oppo- 
site directions from the central engine of a GRB. For an 
observed GRB, the observer should be almost on the jet 
axis. We assume the two jets have the same characteris- 
tics, such as the half opening angle 6'j, the initial Lorentz 
factor Fq, and the kinetic energy Ej/2 (or with equivalent 
isotropic energy Eiso = 2Ej/9'j). Since most afterglow fits 
seem to favor a uniform medium (Panaitescu & Kumar 
2002), we assume a constant particle density n for both 
jets. 

In the above assumptions we have take the standard 
view of a homogeneous jet with sharp edge (Rhoads 1999). 
Recently, a structured jet model (Meszaros et al. 1998; Dai 
& Gou 2001; Rossi et al. 2002; Zhang & Meszaros 2002) 
was proposed, but detailed fits to afterglow data are still 
needed to determine whether it is consistent with current 
observations (Granot & Kumar 2003). Since the standard 
jet model have been well successful in the data fits (e.g., 
Panaitescu & Kumar 2002), we only discuss homogeneous 
jets here. 

2.1. Jet dynamics 

We first discuss the PJ. The PJ evolution can be divided 
into four phases. (1) Initially, F 1/^j, the transverse 
size of the jet is larger than that of the causally connected 
region, therefore the jet evolves as if it was a conical sec- 
tion of a spherical relativistic blast wave. The PJ under- 
goes first a coasting phase with F = Fq. (2) After the jet- 
induced shock sweeps up enough medium, the jet kinetic 

Beijing 100039, China 



2 



Li & Song 



energy is mostly transferred into the shocked medium, and 
the jet begins to decelerate significantly. This occurs at 
the deceleration radius Vd = 3Eiso/ (47rnrQ rUpC^) at an ob- 
server time td = rd/(2rQc), where a relation dr « 2r'^cdt 
has been taken for surperluminal motion. At the decel- 
eration phase, the dynamics could be well described by 
the Blandford & McKee (1976) self-similar solution, where 
r oc r~^/^ oc (3) As the jet continues to decelerate, 

the Lorentz factor drops to l/6j at a time tj = td{ToOj)^^^ , 
corresponding to a radius rj = rd(ro6'j)^/^. Prom the 
time on the transverse size of causally connected regions 
exceeds that of the jet, therefore the sphere approxima- 
tion breaks down, and the jet starts expanding sideways. 
The dynamical evolution in this stage depends on the de- 
gree of sideways expansion. If the lateral velocity in the 
comoving frame equals to the local sound speed, the jet 
will spread quickly with the opening angle increasing and 
F dropping exponentially (Rhoads 1999). At this stage, 
the radius hardly increases and may be regarded as a con- 
stant r « rj, and P oc (4) Finally the jet goes 
into NR phase with P ^ 1 at tpfji = tj/9'j, where the 
evolution can be well described by Sedov- Taylor solution 
/3 oc oc with /? = (1 - l/p2)V2. in summary, 
the evolution is: 
coasting: t <td T = Tq, r oc t 



sphere: 



td<t<tj P oc i"^/^, r oc t^/* 



spreading: tj<t< tjsfR P oc t r ^ rj 

NR: t>tNR /3 oc (r ~ 1), r « t^/s 

Using the above simplified dynamical relation, the NR 
time and radius, defined as where P = 1, are able to be 
calculated as 

1 / 3^;,- ^ / QB'. \ 1/3 
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They are independent of Po and Qj . With typical values 
Ej = W'^Ej^-^x ergs and n = 1 cm~^ (e.g.. Frail et al. 2001; 
Panaitescu & Kumar 2001), we have 
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corresponding radius is rjvij = 2ctNR, therefore for the 
RJ, the NR time is f^jj = Inr + 2rNR/c = ^Inr = 
8.4(£^j-5i/n6'2_^)i/3 yj.. These may be regarded as up- 
per limits, while the values in eqs. (2-4) for the fastest 
spreading cases as lower limits. 

2.2. Non-relativistic phase emission 

For the detailed discussion of afterglow emission in NR 

phase, the readers might refer to Frail, Waxman & Kulka- 
rni (2000), Livio & Waxman (2000) and Waxman (2004), 
we here only give brief introduction. 

Consider also the PJ emission first. After t^R, the evo- 
lution is described by the spherical NR self-similar evolu- 
tion, with the P.J radius given by r = 7'7v_r(V^w_r)^^^i said 
the swept-up particle number is — (4/3)7rnr'^. Sim- 
ply assume that the shocked material forms a thin and 
uniform shell, with the width A = r/r] <^ r and the post- 
shock thermal energy density U = £'j/(47rr^A). The elec- 
tron and magnetic-field energy densities are assumed to 
be constant fraction Cg and cb, respectively, of the post- 
shock thermal energy density U. The electrons are be- 
lieved to be accelerated to form a power law distribution, 
due/d'ye oc for 7e > 7m. Based on the above assump- 
tions, it can be found that the magnetic field strength and 
the minimum electron Lorentz factor evolve as i? oc t~^^^ 
and 7to oc t~^/^, respectively. 

The accelerated electrons will give rise to synchrotron 
emission. The characteristic synchrotron frequency of elec- 
trons with 7j„ is Vm = eB7g/(27rmec), i.e.. 



1/2 

-l^B,-l 



^1/2 f J_ 
\tNR 



GHz 
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where €e,-i = Ce/lO"^, and €3,-1 = es/lO"^. The 
peak specific luminosity is L^^ = NeP/vm where P = 
(4/3)(TTC7m(-B^/87T) is the synchrotron power of an elec- 
tron with 7m. The observed radio frequency is usually 
beyond Vm at NR stage, the specific luminosity is there- 
fore = L^^(i//i/m)(l-P)/2 OC f21/10-3p/2^ ^^^^^ 



Note that this t^R refers to the PJ. 

As for the RJ, it also transits into NR phase at a radius 

TATfl. However, due to the light travel time the observer 
time epoch is delayed by a time 2rNR/c, so the observed 
NR time for the RJ is 

^n'^r = ^NR H = 5tjvfl, (4) 

c 

5 times that of the PJ^. 

There are uncertainties about the lateral spreading of 
jets, and it might be that little spreading occurs, as shown 
in some numerical simulations (e.g., Cannizzo et al. 2004). 
If we assume the extreme case that the lateral velocity is 
zero, then the jet remains conical geometry and continues 
to evolve as if it was a conical section of a spherical rel- 
ativistic blast wave after T < l/9j. In this case, the NR 
time for the PJ is Inr = {l/2){3Ej/27rmpC^ne])^/^ = 
610(i;j,5i/n6l|_i)V3 d, with Oj-i = Oj/lO-'^, and the 

^ It can be shown that the factor 5 is also valid for the wind case where 
source (Dai & Lu 1998; Chevalier & Li 2000). 
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where p — 2 has been taken. The NR light curve flattens, 
compared with the decay, oc t^''', in the spreading phase 
(Rhoads 1999; Sari, Piran & Halpern 1999). 

As for the RJ, when it is still relativistic the emis- 
sion is always beamed away from us and hence invisible. 
At the time t^jj^ after the GRB, the RJ becomes non- 
relativistic and the emission becomes isotropic. The ob- 
served RJ luminosity at this point should be comparable 
to that of the PJ at Inr, because of the same charac- 
teristics of the PJ and RJ. If neglect the NR jet velocity 
with respect to the light speed, the time delay between 
the emission from the PJ and RJ could be roughly fixed 
as 2r/c ~ 2rNR/c = AInr (for r is not rNR)- The RJ 
specific luminosity, then, could be written as, for t > t^'^, 



L^-'{t)c^L^{t-AtNR) cx 



tNR J 

the external density decreases as the square of the distance from the 
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A schematic plot is given in figure 1 for both PJ and RJ 
radio light curves. After rapid increase in flux, the RJ 
emission reaches the peak at around i^'^, and then de- 
creases as a delayed light curve of the NR-phase PJ, but 
showing a steeper decline in the logarithmic plot. Note 
that the transition to NR phase might be gradual, and 
the jet edge might be not sharp, therefore both these can 
smooth the RRB feature on the plot. 

3. OBSERVATION 

3.1. GRB 970508 and GRB 980703 

As discussed above, if a GRB can be monitored for a 
long time in radio bands, e.g., a few years, the RRB may 
be able to be detected. So far, there arc two GRBs that 
have been reported in the literatures with radio monitor- 
ing longer than one year: GRB 970508 and GRB 980703 
(Frail, Waxman & Kulkarni 2000; Berger, Kulkarni & Frail 
2001; Frail et al. 2003). 

Frail, Waxman & Kulkarni (2000) report an extensive 
monitoring of the radio afterglow of GRB 970508, lasting 
450 days after the burst. In the data analysis they found 
that the spectral and temporal radio behavior indicate a 
transition to NR expansion at Inr ~ 100 d, therefore, we 
should expect the RRB at around day 500. However, this 
is unfortunately behind the last observation of this burst, 
so we might have missed that apparent phenomena. 

GRB 980703 was monitored even longer, up to ~ 10^ d, 
as shown in figure 1 of Frail et al. (2003). Two obvious 
flattening appear in the radio light curves: while the late- 
time transition to a constant flux is thought to be the 
host galaxy contribution, the earlier flattening at « 40 d 
is attributed to the transition into NR expansion (Berger, 
Kulkarni & Frail 2001). Modelling the F evolution has 
also inferred a similar value of tNR = 30 — 50 d (Frail et 
al. 2003). Therefore, the expected RRB should arise at 
~ 200 d, and then after another time interval of 200 d, 
around 400 d from the burst, the radio flux should decline 
to a value similar to just before the RRB. At first glance of 
figure 1 of Frail et al. (2003), the radio light curves, some- 
what more extensive at 8.5 and 4.8 GHz, seem to show 
only decays without any re-brightening. However, we no- 
tice that there is no observation data collected between 
210 d and ~ 400 d, therefore the RRB may be missed 
once again due to the lack of observation. 

3.2. Prediction for GRB 030329 

GRB 030329 is the nearest cosmological burst [z = 

0.1685), and has a bright afterglow at all wavelengths. 
Berger et al. (2003) have reported its bright radio after- 
glow up to ~ 70 d after the burst. Since no radio flat- 
tening (due to transition to NR evolution) appears yet^, 
we should calculate Inr using Ej/n derived from available 
data. 

The early breaks (~ 0.5 d) in the R band (Price et 
al. 2003) and X-ray (Tiengo et al. 2004) light curves infer a 

small jet opening angle, leading to a small jet-corrected en- 
ergy release in GRB 030329, which is more than one order 
of magnitude below the average value around which most 

^ Although the X-ray afterglow shows a flattening around day 37 (Tiengo et al. 2004), we believe this is not relevant to the relativistic to 
NR transition of the radio jet, because of no simultaneous radio flattening (Berger et al. 2003), and also conflicted with the still superluminal 
expansion in the radio angular size measurement (Taylor et al. 2004). The X-ray flattening may need other explanations, e.g. an IC component; 
NR transition for the narrow jet in two-component model (Berger et al. 2003). 



GRBs are narrowly clustered (Frail et al. 2001; Panaitescu 
& Kumar 2001). Two models have been proposed to 
solve this problem: a refreshed-shock model with 10 times 
more energy injected later (Granot, Nakar & Piran 2003); 
and a two-component model with most energy in a wider 
jet component (Berger et al. 2003). Both models have 
the same result: the kinetic energy corresponding to the 
late-time radio afterglow is typical for all GRBs. Based 
mainly on the early multi-frequency data, a typical value 
for medium density, n = I cm^'^, was also derived by 
Willingale et al. (2004). Therefore, the estimated value of 
tNR should be typical as in eq. (3), hence tNR ~ 130 d. 
In fact, the radio angular size measurement of this burst 
shows relativistic expanding 83 d after the burst (Taylor 
et al. 2004), consistent with the estimation here. 

In addition, Berger et al. (2003) have inferred Ej^^i/n = 

— 1/2 

O.lAu^ from a snapshot spectrum at tj « 10 d, where 
t'c = fc/10^^ Hz is the cooling frequency. This parameter 
value leads to a too early NR phase, tNR. « 68 d (from eq. 
3), conflicted with the observed superluminal expansion 
(Taylor et al. 2004). However, in the parameter derivation, 
the assumption of jet break time at tj « 10 d might be in- 
correct, because the millimeter observation shows that the 
light curves have already steepened at w 5 d at both 100 
and 250 GHz (Sheth et al. 2003), implying that tj < 5 d. 
Moreover, the angular size evolution can also give con- 
straints to the model parameters. For GRB 030329, Gra- 
not, Ramirez-Ruiz & Loeb (2004) derive Ej^^i/n ~ 0.8 
(sec also Oren, Nakar & Piran 2004) from the observation 
by Taylar et al. (2004). This yields a later transition at 
tNR - 120 d. 

We will take the more reasonable value tNR « 120 d 
(It should be noticed that this might still be the lower 
limit since the most rapid sideways spreading has been 
assumed), and hence the expected RRB should arise at 
~ 1.7 yr after GRB 030329, with the fluxes comparable 
to that at the time Inr- Assuming the same temporal 
slope for late time, an extrapolation of the available radio 
light curves (Berger et al. 2003) to ~ 120 d gives a flux 
~ 0.6 mJy at 15 GHz, or ~ 0.3 mJy at 44 GHz. Dense 
monitoring campaign around ^ 1.7 yr after GRB 030329 
is required to obtain a well-observed RRB profile. 



4. SUMMARY AND DISCUSSION 

In this work we suggest that RRBs are common for 
GRBs with double-sided jets, which come from the RJ 

transition to NR evolution. If we assume the same prop- 
erties for both PJ and RJ and also for the ambient 
medium on both sides, the RRB would arise at a time 
^NR = ^tNR = l-8(i?j,5i/n)^/^ yr after the burst, with a 
flux comparable to that at the time Inr- After tjv/j, the af- 
terglow behaves as a delayed emission of that behind tNR, 
with the time lag of 4:tNR (eq. 7). Unfortunately, no radio 
data could be coflected when RRBs occur for GRB 970508 
and GRB 980703, two longest observed GRBs so far. We 
suppose the RRB of GRB 030329 around 1.7 yr after the 
burst with ~ 0.6 (0.3) mJy at 15 (44) GHz, and urge dense 
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monitoring campaign during that time. It should be no- 
ticed tliat for weak jet-spreading, the RRB might be more 
delayed since t^'^ is much larger (see in the end of §2.1). 

There has been growing evidence for collimated jets in 
GRBs over the past several years, which is coming mainly 
from observations of achromatic breaks in the afterglow 
light curves (e.g., Kulkarni et al. 1999; Stanek et al. 1999). 
However, there are still other explanations for the light 
curve breaks, for example, the transition from the rela- 
tivistic to NR phase of the blast wave at a few days due 
to highly dense medium (Dai & Lu 1999; Wang, Dai & Lu 
2000); the effects of IC scattering flattening or steepening 
the light curves (Wei & Lu 2000); a sudden drop in the 
external density (Kumar & Panaitescu 2000); a break in 
the energy spectrum of radiating electrons (Li & Cheva- 



lier 2001). On the other hand, though black hole-disk 
system with twin jets is generally assumed in GRBs, mag- 
netized, rapidly rotating neutron stars remain contenders, 
and an off-center dipole could lead to a one-sided jet. Since 
RRBs arc only associated with collimated, double-sided, 
relativistic outflows, the detection of RRBs would provide 
straightforward evidence of double-sided jets in GRBs, and 
prove the black hole-disk system formation in the cores 
of progenitors. This make RRB observations significantly 
important. 

Z. Li thanks R. F. Shen for discussions, and L. J. Gou 

for comments. This work was supported by the National 
973 Project and the Special Funds for Major State Basic 
Research Projects. 
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Fk;. 1. — The schematic log- log plot for the radio light curve at u > Pm- The solid and dashed lines correspond to the contribution from 
the PJ and RJ, respectively. The time scalings arc marked, assuming p = 2. After a steep decline in the spreading phase, the PJ light curve 
flatten at tjvjj. The RJ emission rapidly increases to a peak around i^'^ = 5tjvi{, and then declines following the way of the PJ unless a time 
delay of Atpfn, but exhibiting steep light curve in the log-log plot. The dashed-dot line indicates L^ltf^n) = L^'^ (t^'ii)- 



